Assessment of Tunisian Trichoderma Isolates on Wheat Seed Germination, Seedling Growth and Fusarium Seedling Blight Suppression

Beneficial microorganisms, including members of the Trichoderma genus, are known for their ability to promote plant growth and disease resistance, as well as being alternatives to synthetic inputs in agriculture. In this study, 111 Trichoderma strains were isolated from the rhizospheric soil of Florence Aurore, an ancient wheat variety that was cultivated in an organic farming system in Tunisia. A preliminary ITS analysis allowed us to cluster these 111 isolates into three main groups, T. harzianum (74 isolates), T. lixii (16 isolates) and T. sp. (21 isolates), represented by six different species. Their multi-locus analysis (tef1, translation elongation factor 1; rpb2, RNA polymerase B) identified three T. afroharzianum, one T. lixii, one T. atrobrunneum and one T. lentinulae species. These six new strains were selected to determine their suitability as plant growth promoters (PGP) and biocontrol agents (BCA) against Fusarium seedling blight disease (FSB) in wheat caused by Fusarium culmorum. All of the strains exhibited PGP abilities correlated to ammonia and indole-like compound production. In terms of biocontrol activity, all of the strains inhibited the development of F. culmorum in vitro, which is linked to the production of lytic enzymes, as well as diffusible and volatile organic compounds. An in planta assay was carried out on the seeds of a Tunisian modern wheat variety (Khiar) by coating them with Trichoderma. A significant increase in biomass was observed, which is associated with increased chlorophyll and nitrogen. An FSB bioprotective effect was confirmed for all strains (with Th01 being the most effective) by suppressing morbid symptoms in germinated seeds and seedlings, as well as by limiting F. culmorum aggressiveness on overall plant growth. Plant transcriptome analysis revealed that the isolates triggered several SA- and JA-dependent defense-encoding genes involved in F. culmorum resistance in the roots and leaves of three-week-old seedlings. This finding makes these strains very promising in promoting growth and controlling FSB disease in modern wheat varieties.


Introduction
Organic and chemical-free farming methods to control biotic pathogens and enhance yield are increasing in popularity and practice currently due to the growing concern about Microorganisms 2023, 11 the potential negative impact of synthetic chemicals used in conventional farming. This improved yield has been accompanied by modern breeding practices for half of the 20th century in order to select the most effective crop genotypes that take full advantage of farming inputs, thus limiting any potential added value coming from beneficial plantmicrobe interactions [1,2]. Due to this, there is an urgent need to develop promising new ways to increase both crop performance and tolerance to biotic and abiotic factors without harming human health or the environment. Worldwide, wheat is one of the major crops, providing food for about 35% of the world's population and making up 17% of the global cultivated area [3]. In Tunisia, cereals and their derivatives constitute the base of all diets and the main source of calories. Tunisia's average rate of wheat consumption per capita is 184 kg per year, making Tunisia one of the largest wheat consumers in the world, a rate which is expected to remain at the same level in the coming years [4,5].
Cereals are susceptible to fungal diseases such as Fusarium seedling blight (FSB) caused by Fusarium fungi, including Fusarium culmorum, the dominant pathogen in both Tunisia and the Mediterranean region [6]. FSB is responsible for extensive damage to growing seedlings; it can cause significant yield losses and impact the quality of the harvested grain resulting in reduced market value [7]. This disease is most severe in cool and wet conditions, which are common in the winter in Tunisia [8]. Symptoms of FSB include stunted seedlings, yellow leaves and root rot, as well as infected seeds which fail to germinate or die soon after germination. In some severe cases, the entire seedling can collapse and rot [7,9]. F. culmorum is a soil-borne fungus that can persist in the soil as chlamydospores for many years, and it can also survive in stubble residues of cereals and other grasses as hyphae [8,10]. This ability makes F. culmorum difficult to manage and has major implications for designing effective strategies for disease management. F. culmorum is capable of causing both FSB and Fusarium Head Blight (FHB) in wheat. FHB disease occurs in the spikes and is characterized by visible symptoms, such as reddish and scabby spikes, which can result in the accumulation of mycotoxins in the grains [11,12]. Control of FSB has received less attention than control of Fusarium head blight (FHB) because of the risk of mycotoxin contaminations of grain, even though FSB could provide a pathogen source for subsequent epidemics of FHB [12,13]. Control of FSB could be achieved by chemical, cultural and biological strategies [14,15]. Using pesticides as a seed treatment could be problematic because of the toxic residues associated with fungicides [16,17]. Some cultural practices may help reduce the fungal population in the soil, such as deep plowing or removing crop residues, but the effectiveness of these measures is limited due to the persistence of the pathogen in the soil [18]. Another management strategy, crop rotation, is not generally practiced in Tunisia because the market has always had a significant influence on the diversity of crops grown in the country. The aim of increasing food self-sufficiency has led to the promotion of certain crops; durum wheat, for example, is a dominant crop taking approximately 60% of the cereal growing area [19]. In addition, control of FSB has been challenging both due to the lack of fully resistant wheat cultivars and because Tunisia is now considered a hot spot for climate change [20].
The use of beneficial microorganisms which are isolated from soil or plant tissues to control plant pathogens and promote plant growth seems to be a very promising and cost-effective component of an integrated disease management plan for sustainable modern agriculture [21].
From this perspective, Trichoderma-based products are a segment that is growing in popularity in the agriculture industry, as Trichoderma spp. can help improve soil health and plant growth. These products are used for a variety of applications, including soil inoculants, seed treatments, biopesticides and biofertilizers. Indeed, Trichoderma spp. are common rhizospheric inhabitants, and the ability of certain species to colonize plant roots allows them to engage with their host in complex molecular dialogues that promote plant growth and productivity [22]. In addition, Trichoderma spp. are considered BioControl Agents (BCAs), conferring protection against a wealth of plant pathogens [23]. Multiple molecular mechanisms have been suggested to be responsible for its biocontrol abilities, such as interference with the life cycles of plant pathogens by mycoparasitism (or hyperparasitism); using the secretion of antibiotics (volatile organic compounds, VOCs) and extracellular cell wall-degrading enzymes (CDWEs); competition for nutrients and space; chemical modification of the environmental conditions; and modulation of the host's innate immunity performance [24][25][26][27][28]. The local and systemic elicitations of the plant defensive responses are correlated with a significant increase in disease resistance against a broad spectrum of plant pathogenic fungi, such as various Fusarium spp. [29][30][31][32][33][34].
Some Trichoderma strains were isolated from Tunisian soil and plant habitats and were studied as beneficial agents in wheat [33,35]. However, the diversity of available beneficial strains is very limited, predominantly focused on some species that primarily belong to one of the five taxonomic sections of the "Trichoderma harzianum species complex" (THSC) [36][37][38]. This gap concerns all of the plant species of agronomic interest. In this respect, exploring the greater diversity of Trichoderma and evaluating their biocontrol abilities would be of great benefit in minimizing the pathogens spread in wheat. However, the success of a fungal biocontrol agent depends not only on the complex interactions it is able to establish with the plant host and its related pathogens but also on its aptitude to colonize a rhizospheric environment modeled by the crops already present, the soil structure and the remarkable endogenous microbial diversity.
Based on the hypothesis that ancient wheat genotypes could represent a source of beneficial fungi such as Trichoderma, we focused our research on Florence Aurore, an ancient wheat variety that was cultivated in an organic farming system in Tunisia, for this study. Our aim was to check this variety for the presence of beneficial Trichoderma species in the rhizosphere through isolation, molecular characterization (ITS, tef1 and rpb2) and assessment of their plant growth promotion abilities (with the production of indolelike compounds, hydrogen cyanide, ammonia and lytic enzymes, or the solubilization of phosphate). In addition, the isolated strains were evaluated for their possible biocontrol activity against F. culmorum, in the context of Fusarium seedling blight, through dual confrontation assays and assessments of diffusible and volatile compound productions. Finally, to confirm these PGP and biocontrolling features, in planta experiments were conducted in greenhouse conditions on Khiar, one of the most marketed wheat varieties in Tunisia, followed by the quantification of defense-related gene expression in the leaves and roots to assess the potential of the Trichoderma isolates as plant immunity inducers. Our data provide a set of new Tunisian Trichoderma strains that could contribute to the use of field applications of Trichoderma fungus as a biocontrol agent.

Tunisian Sampling Location
Field sampling was conducted in May 2020. The study site was located in "Gousset el bey" in the governorate of Bizerte at 36.924096 • N, 9.699089 • E, in the northwestern part of Tunisia ( Figure S1 in Supplementary Materials). This area is characterized by a warm Mediterranean climate with hot, dry summers and humid and cool winters. The average annual temperature range is 18.4 • C, and the rainfall is around 547 mm per year (https: //en.climate-data.org/africa/tunisia/bizerte/bizerte-3551/; accessed on 31 May 2023).

Isolation and Culture Conditions
Trichoderma strains were isolated from rhizospheric soil of healthy plants of the wheat variety Florence Aurore cultivated in an organic farming system in Tunisia. Rhizospheric soil samples were randomly collected at grain-filling stage, i.e., Growth Stage 87, according to the code defined by Zadoks et al. [39]. One gram of soil tightly adhering to the roots was incubated in 10 mL of sterile Deshydrate Bouillon Potato Dextrose (Potato Dextrose Broth, PDB) in a test tube and shaken at 200 rpm for 30 min at 28 ± 2 • C. From this starting solution, serial dilutions were prepared at the following concentrations, 10 −1 , 10 −2 and 10 −3 , and then plated on Trichoderma-Selective Medium (TSM) [40]. After 5 days of The phylogenetic analyses were conducted on the concatenated gene version 5 -ITS (575 bp)/tef1 (1260 bp)/rpb2 (930 bp)-3 , with the reference sequences downloaded from GenBank that correspond to the most related sequences. All of the nucleotide data were weighted equally, and gaps were treated as missing characters. After removing any ambiguously aligned regions, the alignments were 2639 nucleotides. The multiple DNA sequence alignments were performed with the Clustal-W algorithm. The maximum likelihood analysis was employed to construct a phylogenetic tree as implemented in the PhyML program. The topology of the tree was evaluated by bootstrap analysis of the sequence data based on 1000 random resamplings.

In Vitro Confrontation Assay
The antagonism of Trichoderma strains against F. culmorum was evaluated with dual culture assays, according to Skidmore and Dickinson [46]. Mycelial plugs of 5 mm diameter of F. Culmorum and Trichoderma strains (7-day-old) were placed on opposite sides of PDA Petri plate (8.5 cm diameter, containing 20 mL of PDA medium) at 5 cm from each other and at equal distance from the periphery. Co-cultures were incubated at 28 • C for six days. The radial growth of the F. culmorum in the presence of Trichoderma strains was measured, and the inhibition rate (IR) was calculated as follows: IR = R1 − R2/R1 * 100, where R1 is the colony radius of F. culmorum in the control plate, and R2 is the colony radius of F. culmorum in the presence of Trichoderma strain. Each biological condition was repeated five times.

Determination of Diffusible and Volatile Metabolite Productions
The activity of Trichoderma diffusible compounds was tested using cellophane method described by Dennis and Webster [47]. Strains were grown on Petri dishes containing a sterile cellophane sheet over the PDA medium. After five days of incubation at 28 • C, cellophane was removed, and a 5 mm plug of F. culmorum was placed in the center of each Petri plate. F. culmorum grown on untreated PDA served as control. Assays were carried out 5 times, and plates were incubated at 28 • C for five days. Results were expressed as the inhibition rate (IR) calculated according to the formula mentioned above, where R1 is the colony radius of F. culmorum in the control plates, and R2 is the colony radius of F. culmorum in the treated plates. To verify the volatile metabolite production, a qualitative method was used. Before inoculation, a 1.5 cm wide agar strip was removed from the mid-portion of PDA medium. A 5 mm diameter disk of Trichoderma strain was placed on one side, while a disk of F. culmorum was placed on the opposite side. PDA plates inoculated only with F. culmorum served as controls. All plates were sealed with parafilm, and this assay includes 5 replicates per biological condition. After seven days of incubation at 28 • C, the RI was determined according to the abovementioned formula.

Qualitative Determination of Detoxifying and Lytic Enzyme Production
Qualitative determination of extracellular enzymes was performed on 7-day-old Trichoderma cultures. All assays were carried out in triplicate, and a non-inoculated plate was used for each biochemical test as negative control.
Catalase. Catalase activity was assessed with 3% v/v hydrogen peroxide. A 5 mm diameter plug of actively growing Trichoderma culture was put on a glass slide in a Petri plate, on which a drop of hydrogen peroxide was applied. The appearance of immediate bubbling reveals the positive catalase ability of the strain.
Protease. For protease enzyme detection, Trichoderma strains were grown on Glucose-Yeast-Peptone medium (GYP) with the following composition (per liter of distilled water): glucose, 1 g; yeast extract, 0.1 g; peptone, 0.5 g; agar, 15 g, supplemented with 1% Skim milk (pH 6.5) [48]. After seven days of incubation at 28 • C, a positive result for protease production was revealed by the appearance of a clear zone around the Trichoderma colony.
Amylase. To evaluate the amylase activity, Trichoderma strains were inoculated following the previous conditions on GYP medium enriched with 1% soluble starch as substrate for amylase enzyme [49]. After seven days of incubation, plates were flooded with an aqueous solution of Iodine 1% (w/v) in potassium iodide 2% (w/v) for 10 min, and then, the iodine solution was carefully removed by decanting each plate. Amylase activity was revealed by the presence of a clear halo around Trichoderma colony.
Chitinase. Chitinase assay was performed by growing Trichoderma strains on a basal medium comprising (per liter of distilled water): colloidal chitin, 1.5 g; K 2 HPO 4 , 2.7 g; MgSO 4 ·7H 2 O, 0.7 g; NaCl, 0.5 g; KCL, 0.5 g; yeast extract, 0.13 g, agar 15 g, and pH 5.5 [50]. After seven days of incubation, according to the previously cited conditions, plates showing a transparent halo around the colony were considered as Chitinase positive. The method of Roberts and Selitrennikoff [51] was used for colloidal chitin preparation. Five grams of chitin powder was added to 60 mL concentrated HCl and kept at 4 • C overnight with vigorous stirring. The mix was added to 500 mL of ice-cold 95% ethanol with rapid stirring and left overnight at room temperature. The precipitate was collected by centrifugation at 5000× g for 20 min at 4 • C, washed with sterile distilled water several times until the colloidal chitin became neutral (pH 7.0) and stored in the dark at 4 • C before being used. Endoglucanase.
Czapek-agar medium was used to reveal cellulose degradation [52]. It is composed of (per liter of distilled water) NaNO 3 , 2 g; KH 2 PO 4 , 1 g; KCl, 0.5 g; MgSO 4 ·7H 2 O, 0.5 g; FeSO 4 ·7H 2 O, 0.01 g; agar, 20 g, and supplemented with 0.1% carboxymethylcellulose (CMC) as the sole carbon source (pH 4.5 adjusted with 100% glacial acetic acid). Petri dishes were inoculated with plugs of the different Trichoderma strains and incubated at the same cited conditions. After seven days, plates were flooded with an aqueous solution of Congo red (0.1% w/v). After five minutes of reaction, this solution was gently removed, and dishes were washed with 5 M NaCl to reveal the degradation halo. The capacity to solubilize phosphate was tested on solid Pikosvskaya medium (PVK) [53] that had the following composition (per liter): 10 g of glucose, 0.5 g of yeast extract, 0.5 g of (NH 4 ) 2 SO 4 , 0.1 g of MgSO 4 , 0.2 g of KCL, 15 g of agar and 5 g of Ca 3 (PO 4 ) 2 as sole source of phosphate. Plates were inoculated with a plug of 5 mm of Trichoderma (7 days old) with 3 repetitions for each strain per plate. After incubation for 7 days at 28 • C, the presence of a clear halo around the colony indicates the presence of phosphate dissolving activity. The experiment was conducted twice, with three sets of replication plates.

Colorimetric Detection of Indole-Related Compounds
Indole-related compound production was determined using a colorimetric method [54]. Two plugs of 5 mm of each Trichoderma strain (7 days old) were transferred in 10 mL of Potato Dextrose Broth (PDB) supplemented with 0.1% L-Tryptophan and then incubated under shaking at 28 • C for 72 h. Supernatant was collected from 1 mL of culture after 5 min of centrifugation at 14,000 rpm, then 100 µL of this supernatant was added to 200 µL of Salkowski reagent in triplicate. After incubation at room temperature for 30 min, the optical density was recorded at 530 nm [55]. The amount of indole-related compounds was calculated from a standard curve generated by suspending IAA in 100% ethanol at a concentration of 1 mg/mL and then diluted in PDB medium to a concentration of 1, 2, 3, 4, 5, 7, 8, 9 and 10 µg/mL. PDB with or without L-Tryptophan served as controls. The experiment was conducted with three biological replicates.

Hydrogen Cyanide (HCN) Production
Trichoderma spp. was grown on solidified PDA supplemented with 4.4 g/L of Glycine or Succinate. White filter paper discs cut into the same size as the upper lid of the Petri dish were immersed in an alkaline Picric acid solution (Picric acid, 0.5% (v/v), Sodium carbonate, 1.25% (w/v); in water), and carefully positioned on the lid of each plate. Plates were sealed with Parafilm and incubated for seven days at 28 ± 2 • C. After incubation, HCN production was observed by the color changes of the filter paper from yellow to light brown or reddish brown, which indicated the production of HCN [56]. Because HCN can be produced at very low concentrations, the filter paper was then immersed in 5 mL of distilled water to dissolve the precipitate, and the optimum density was measured at 625 nm using a spectrophotometer [57]. Each measurement corresponds to the mean of five biological replications.

Ammonia Production
Ammonia production was determined using a colorimetric method [58]. Each Trichoderma strain was tested for the production of ammonia in peptone water. Two plugs of 5 mm of actively growing Trichoderma (7 days old) were added to 10 mL peptone water and then incubated under shaking at 28 • C for 7 days. Cultures were centrifuged at 3500 rpm for 10 min. A volume of 50 µL of Nessler's reagent (7% KI; 10% HgCI 2 ; 50% aqueous solution of NaOH (32%)) was added to 1 mL of every culture supernatant. Development of yellow to brown color indicates the presence of ammonia. The optical density was measured at 450 nm, and ammonia production was calculated using a standard curve of ammonium sulfate ((NH 4 ) 2 SO 4 ) solution at different concentrations of 0.1, 0.5, 1, 1.5, 2, 2.5, 3, 3.5 and 4 µmol/mL. Each measurement corresponds to the mean of three biological replications.

Bioassay of Trichoderma Strains for Plant Growth Promoting and Fusarium Seedling Blight Disease Suppression
The plant growth-promoting effect, as well as the suppression of Fusarium seedling blight (FSB) disease of 6 isolated Trichoderma strains on durum wheat, was tested by germinating the seeds in axenic conditions and then by monitoring their growth under greenhouse conditions for 3 weeks. The variety "Khiar", which is one of the most used and productive durum wheat varieties in Tunisia, has been chosen for this study [59]. Seeds were provided by the National Institute of Agronomic Research of Tunisia (INRAT).

Effect on Seed Germination in Axenic Conditions
Seeds were surface sterilized as described by Fernandez and Chen (2005) [60] with minor modifications. Seeds were soaked in an aqueous solution of 0.5% sodium hypochlorite (NaClO) for 5 min, washed three times with sterile distilled water, dried overnight on filter paper under aseptic conditions and placed in sterile dishes for further use. For fungal preparation, Trichoderma strains were grown on PDA medium under 12 h regime of fluorescent light at 26 ± 2 • C for 7 days to promote sporulation. Fungal spore suspension of F. culmorum was prepared from a 2-week-old culture. The collected fungal suspension was filtered, macrospores were counted, and the concentration was adjusted to 10 7 spores/mL in sterile PDB.
Seeds were soaked for 24 h under slow shaking in 20 mL of Trichoderma-PDB suspensions (10 7 spores/mL) to promote Trichoderma growth on seeds (50 seeds in a falcon tube of 50 mL for each repetition). Seeds soaked only in PDB served as control, and each treatment was repeated 4 times. Seeds were dried on filter paper under aseptic conditions and then sown in square Petri dishes of 12 * 12 cm (50 seeds per plate) containing 1% of water agar medium (10 g/L) previously autoclaved and sterilized for 30 min under UV radiations [61]. For Fusarium seedling blight suppression, seeds were prepared in a similar way, coated with Trichoderma strains and then sprayed with 1 mL of F. culmorum spore suspension adjusted at 10 7 spores/mL. Seeds sprayed only with F. culmorum served as positive control. All the following treatments, i.e., for non-treated seeds, seeds treated only with Trichoderma, Microorganisms 2023, 11, 1512 8 of 36 seeds treated with Trichoderma and F. culmorum and seeds treated only with F. culmorum, were repeated 4 times. All plates were sealed with parafilm and incubated in the dark at 25 • C for 5 days, and some parameters, which include the final germination percentage (FGP), fresh weight of seedlings and root and shoot lengths, were measured. FGP was calculated using the following formula: FGP = germinated seeds/total seeds * 100.

Effect on Seedling Growth under Greenhouse Conditions
After 7 days of germination in axenic conditions, PDB control and fungal-treated seedlings (treatments described in the Section 2.5.1) were carefully removed from water agar medium and transferred into pots of 6 * 6 * 6 cm, with 1 seedling per pot. Each experiment incorporates six biological repetitions. Substrate consisted of 1/3 sterilized sand and 2/3 potting soil. All pots were placed in greenhouse under the following conditions: 16/8 h light/dark photoperiod under an illumination of 150 µmol/m 2 /s photosynthetic photon flux density, a temperature of 23 ± 0.5 • C, and a relative humidity of 65 ± 5%. Pots were watered every two days with tap water for 21 days until all plants were harvested to complete all the ecophysiological measurements (fresh and dry weights, shoot and root lengths, cf Section 2.5.3) and molecular analyses (cf Section 2.7).

Measurements of Total Chlorophyll, Epidermal Flavonols, Nitrogen Balance Index (NBI)
Chlorophyll content (Chl), epidermal flavonols (Flav) and NBI were measured in vivo using Dualex sensor (Force-A, France), which is a clamp that performs instantaneous and non-destructive quantifications [62]. All proxy measurements were done on uniform, fully developed and light-exposed leaves of the same plant 21 days after fungal treatment.
2.6.4. Total RNA Isolation, cDNA Synthesis, and Quantification of the Defense-Encoding Gene Accumulation From the experimental set-up used for the ecophysiological analyses and concomitantly to the execution of the ecophysiological measurements, the aerial and root tissues were collected from fifteen randomly selected plants, which were then distributed into 3 biological replicates, each comprising the tissues of 5 plants. All samples were immediately frozen in liquid nitrogen and stored at −80 • C before being used. Samples were ground to a fine powder with 3 metallic balls for 90 s (Ø 5 mm) in Eppendorf tubes. Total RNA was extracted using a CTAB extraction buffer (Cetyltrimethylammonium Bromide), including the elimination of the genomic DNA as previously described in Ben Amira et al. (2018) [63]. The quality and quantity of the recovered total RNA were estimated by a NanoDrop ND-1000 spectrophotometer. Two µg of total RNA were reverse-transcribed with Oligo-dT using the SuperScript ® III First-Strand Synthesis System for RT-PCR (Invitrogen, Carlsbad, CA, USA) according to the manufacturer's instructions. The quantification of the transcript levels of genes was measured by real-time quantitative PCR from 1 µL of cDNA diluted 30-fold with sterile water. Expression analysis of the targeted genes involved in the defense responses (listed in Table S1 in Supplementary Materials) was performed with an Applied Biosystems StepOnePlusTM Real-Time PCR System (Applied Biosystems, Foster City, CA, USA) using TakyonTM Rox SYBR ® MasteMix dTTP Blue (Eurogentec, Liege Science Park, Belgium), according to manufacturer's recommendations. PCR cycling conditions were 1 cycle for 3 min at 95 • C, followed by 40 cycles for 10 s at 95 • C, 15 s at 60 • C, and 10 s at 72 • C. PCR reactions were ended by generating a dissociation curve of 60-95 • C with an increment of 0.3 • C/15 s to ensure primer dimers and nonspecific amplifications. Three reference genes (Actin (ACT), Heterogeneous nuclear ribonucleoprotein (RPN), and RNaseL inhibitor-like (RLI)) were used to normalize the expression results of the target genes. The reference genes were chosen from a panel of widely used reference genes and specifically selected for their stable expressions between samples and treatments. Moreover, each referrer belongs to protein families involved in different cellular processes in order to minimize the risk of co-regulations. The expression level of the normalized target genes in each RNA preparation was calculated using the 2 −∆∆Ct [64]. All reactions were performed in triplicate.
Primer pairs are detailed in Table S1 in Supplementary Materials. PCR efficiencies (E) and conditions were determined by comparing threshold values in a dilution series of the RT product (×10, ×20, ×40, ×80), followed by non-template control for each primer pair. PCR efficiencies are given in Table S1 in Supplementary Materials.

PCR Detection of Trichoderma and F. culmorum Strains in Wheat Roots and Leaves
The presence of Trichoderma or F. culmorum isolates was qualitatively characterized in the cDNA samples used for gene expression defense analysis. Amplifications targeted the tubulin-encoding gene with primer pairs specific to each fungal genus (Trichoderma: 5 -GTACTAAGTTGTTTCTTTGCTGTTG/5 -CTCTTGTACATACACCAATTGCTC; Fusarium: 5 -GAGTACTAAGCGGTTTCGGATGC/5 -GCATGCCGCGTGCGCGCAAGGTTC; (5 -Forward/5 -Reverse)) [63]. The PCRs were carried out on cDNAs diluted by 1:5 and 1:20, according to the PCR cycling conditions described above, and with a hybridization temperature fixed at 68 • C.

PCR Statistical Analysis
Ecophysiological response variables were collected from at least three biological replicates. The standard one-way Analysis of Variance (ANOVA) was carried out to test the effect of Strain factor on each of these variables. Homogeneity of variance between strains was examined using Levene's test [65] using the car::leveneTest function [66], and normality was examined through visual inspection of the residuals as well as the Shapiro-Wilk test [67] using the stats:shapiro.test function using R Statistical Software (v4.2.2; R) [68]. If the homogeneity of variance assumption was violated, Welch One-Way ANOVA [69] was carried out using the stats:oneway.test function. Multiple comparisons of means were carried out using Tukey's HSD post hoc test, with the "rstatix::tukey_hsd()" for standard ANOVA or the Games-Howell test [70] with the "rstatix::games_howell_test()" for Welch ANOVA [71]. Omega squared was used as a measure of the effect size for standard ANOVA and as estimation for Welch ANOVA [72]. It is an estimate of how much variance in the response variables is explained by the Strain factor. For omega squared, 0.01 to <0.06 indicates a small effect size, 0.06 to <0.14 is a medium effect size, and ≥0.14 is considered a large effect size [73]. It was computed with MOTE::omega. F function [74]. The Kruskal-Wallis test [75] based on ranks of data was computed with stats:kruskal.test function and used instead of a one-way ANOVA in case the normality assumption did not hold. The rank epsilon squared was then computed with rcompanion: epsilonSquared and utilized as a measure of the Kruskal-Wallis test effect size [76]. It varied between 0 and 1 and indicated the ratio of variance in the response variable explained by the explanatory variable, the Strain factor. It also indicated the degree to which one strain had data with higher ranks than the other strain. For rank epsilon-squared, 0.01 to <0.08 indicates a small effect size, 0.08 to <0.26 is a medium effect size and ≥0.26 is considered a large effect size [77]. Kruskal-Wallis test was followed by Dunn's post hoc test [78] using agricolae::kruskal [79] for multiple comparisons of medians. All statistical analyses were performed using R Statistical Software (v4.2.2; R) [68].

Phylogenetic Position of the Trichoderma Isolates
A total of 111 strains were isolated on TSM medium from rhizospheric soil samples. In order to identify the potential strains from the diversity of colonies, a fragment of approximately 575 bp of the ITS region was obtained from PCR amplification using ITS1 and ITS4 primers, followed by forward and reverse sequencing. The BLASTn alignment on the NCBI database enabled us to retrieve six different sequences that share about 98% to 100% similarity with Trichoderma strains. These sequences were divided into three large groups with highly significant bootstrap values (NJ (Neighbor-Joining) method) ( Figure 1): a first group of 66.66% from the total isolates represented by three different Trichoderma harzianum strains, a second group of 14.41% represented by one Trichoderma lixii and a third group of 18.91% represented by two Trichoderma sp. close to the T. harzianum complex. Every sequence was recorded in NCBI, and the related GenBank accession numbers are given in Figure 1.
[83] packages. Differential gene expressions were selected after false discovery rate control using Benjamini-Hochberg procedure [84] and filtered at 5% error level.

Phylogenetic Position of the Trichoderma Isolates
A total of 111 strains were isolated on TSM medium from rhizospheric soil samples. In order to identify the potential strains from the diversity of colonies, a fragment of approximately 575 bp of the ITS region was obtained from PCR amplification using ITS1 and ITS4 primers, followed by forward and reverse sequencing. The BLASTn alignment on the NCBI database enabled us to retrieve six different sequences that share about 98% to 100% similarity with Trichoderma strains. These sequences were divided into three large groups with highly significant bootstrap values (NJ (Neighbor-Joining) method) ( Figure 1): a first group of 66.66% from the total isolates represented by three different Trichoderma harzianum strains, a second group of 14.41% represented by one Trichoderma lixii and a third group of 18.91% represented by two Trichoderma sp. close to the T. harzianum complex. Every sequence was recorded in NCBI, and the related GenBank accession numbers are given in Figure 1. Interestingly, all six isolates were placed in the T. harzianum complex [36][37][38]. Since ITS loci analyses are insufficient to distinguish the species' delimitation within a closely Interestingly, all six isolates were placed in the T. harzianum complex [36][37][38]. Since ITS loci analyses are insufficient to distinguish the species' delimitation within a closely related Trichoderma species complex such as T. Harzianum [85][86][87], partial regions of loci encoding tef1 and rpb2 were sequenced from these six "T. harzianum" isolates. These loci are powerful sources of information due to their inherent interspecific variations and are extensively used to accurately resolve the taxonomy of Trichoderma [88,89]. The analysis of concatenated ITS-tef1-rpb2 sequences confirmed that the six isolates belong to a common complex. Moreover, it facilitated the identification of the two Trichoderma sp. isolates, as being the T. afroharzianum (Tahz03) and T. lentinulae (Tlen01) species. Similarly, the species' range was substantially clarified for the three T. harzianum, which are precisely affiliated with two T. afroharzianum (Tahz01 and Tahz02) and one T. atrobrunneum (Tatr01). T. atrobrunneum is a new species clade closely related to the T. Afroharzianum taxon but clearly distinguished within the T. harzianum sub-complex [90]. Finally, this multi-loci analysis confirmed the correct identification of T. lixii (Tlix01), which is phylogenetically close to T. lentinulae and forms a highly supported clade with it. These two clades, T. afroharzianum/T. atrobrunneum and T. lixii/T. lentinulae, have been separated from the T. harzianum complex [38,91].
Each strain exhibited the expected morphological characteristics (macro-and microscopic) of the Trichoderma genus ( Figure S2 in Supplementary Materials). On the Potato Dextrose Agar (PDA) medium, Trichoderma strains appeared as fast-growing colonies, concentric halos and floccose or compact surfaces that looked similar to tufts. The mycelium started out white, then turned green due to the production of conidia. Under an optical microscope, abundant sporulation of conidia with smooth or rough appearance was observed, conidiophores were branched and irregularly verticillated, and phialides were generally ampliform or fusiform and arranged in clusters, according to Gam et al. [92].
All of these species have also been reported all across the world; however, to our knowledge, this is the first report of T. atrobrunneum, T. lixii and T. lentinulae in an agroecological context in Tunisia, and they have never been introduced as beneficial agents in wheat.
Some Trichoderma strains were isolated from Tunisian soil and plant habitats. Their characterization was based on evaluating their occurrence and diversity within contrasting bioclimatic zones [93,94] or on their PGP and/or BCA abilities on perennial (T. harzianum on Olive trees [32]) or annual crops (e.g., T. harzianum and T. viride on Potatoes [95] and Faba bean [96]; and Trichoderma spp. on wheat [33,34]). The strains mentioned in these studies represented different ranges of Trichoderma complexes, and these strains and some of our isolates presented common strain complexes, of which T. (afro)harzianum appeared to be the most representative. However, little is known about the diversity of Trichoderma spp. that occurs or interacts with wheat, with most studies focusing on T. harzianum and more marginally on T. gamsii, T. viride, T. koningii or sp. [16,33,35,97,98]. T. harzianum T22, originally derived from the fusion of two auxotrophic strains, T-95 and T-12, and the first fungus registered by the Environmental Protection Agency (EPA) for biocontrol of plant diseases in 1989 [99], is still one of the most widely used biopesticides today [16,100,101]. The genus Trichoderma has more than 150,000 members, and of the hundreds of Trichoderma sp. species recorded, only a fraction is exploited as biocontrol agents on field crops. The diversity of the above-mentioned Tunisian strains added to our own (with the originality of the T. atrobrunneum, T. lixii and T. lentinulae) and remained, nevertheless, a very promising source of beneficial agents. To this end, the study of this fungal genus is still a very active research area, especially as the beneficial influences of these fungi are dependent on the fungal strains and the plant species with which they interact. Likewise, this situation is paradoxically prevalent in durum wheat, especially during FSB on durum wheat caused by F. culmorum. Because the plant growth-promoting potential varies between different Trichoderma species and strains, and also because local isolates are now favored over Trichoderma's subsequent use as biofertilizer and biocontrol agents, we have evaluated the plant growth-promoting abilities and the antagonism activities of these six native Trichoderma strains isolated from wheat rhizosphere against that of F. culmorum in vitro and in vivo.

Trichoderma Isolates Exhibit Differential PGP Abilities
During germination in axenic condition, coating seeds with all six Trichoderma strains showed no significant effect on the final germination percentage effect (p-value for Kruskal-Wallis rank sum test = 0.4971; rank epsilon2 = 0.199) or on total fresh weight of 7-day-old seedlings (p-value for standard ANOVA = 0.246; omega squared = 0.029) (Figure 2a).
In addition, no particular symptom related to any pathogenicity expression of our tested Trichoderma strains was observed on the variety of wheat Khiar used in this study, whether during germination or seedling growth experiments (Figure 2c; Figures S3 and S4 in Supplementary Materials). These first symptomatological inspections are of the utmost significance and should be performed on all wheat genotypes commonly grown in Tunisia. Indeed, some studies have shown that the treatment with Trichoderma either had no effect or was even detrimental. In this respect, several Trichoderma species are included among pathogens of cultivated plants: e.g., T. viride is the causal agent of onion green mold rot [102], and T. afroharzianum is responsible for the occurrence of a new ear rot disease in maize in Europe [103].

Trichoderma Isolates Exhibit Differential PGP Abilities
During germination in axenic condition, coating seeds with all six Trichoderma strains showed no significant effect on the final germination percentage effect (p-value for Kruskal-Wallis rank sum test = 0.4971; rank epsilon2 = 0.199) or on total fresh weight of 7-dayold seedlings (p-value for standard ANOVA = 0.246; omega squared = 0.029) (Figure 2a). In addition, no particular symptom related to any pathogenicity expression of our tested Trichoderma strains was observed on the variety of wheat Khiar used in this study, whether during germination or seedling growth experiments (Figure 2c; Figures S3 and S4 in Supplementary Materials). These first symptomatological inspections are of the utmost significance and should be performed on all wheat genotypes commonly grown in Tunisia. Indeed, some studies have shown that the treatment with Trichoderma either had no effect or was even detrimental. In this respect, several Trichoderma species are included among pathogens of cultivated plants: e.g., T. viride is the causal agent of onion green mold rot [102], and T. afroharzianum is responsible for the occurrence of a new ear rot disease in maize in Europe [103].
Seven days after planting seedlings in pots, the plant growth-promoting ability was clear, and every isolate exhibited a significant (p-value for Welch's ANOVA = 1.53 × 10 −7 ; estimated omega squared = 0.562) biostimulatory effect on shoot height (Figure 2b,c). This impact was still significant 14 days later (p-value for Welch's ANOVA = 1.80 × 10 −9 ; estimated omega squared = 0.693) (except for Tlen01) and even up to 21 days later (p-value for standard ANOVA = 1.75 × 10 −5 ; omega squared = 0.193) (except for Tlen01 and Tatr01). The effect size for the shoot height after 7 days, as well as after 14 days and 21 days, was large; the observed differences among strains for the shoot height were of practical significance. Seven days after planting seedlings in pots, the plant growth-promoting ability was clear, and every isolate exhibited a significant (p-value for Welch's ANOVA = 1.53 × 10 −7 ; estimated omega squared = 0.562) biostimulatory effect on shoot height (Figure 2b,c). This impact was still significant 14 days later (p-value for Welch's ANOVA = 1.80 × 10 −9 ; estimated omega squared = 0.693) (except for Tlen01) and even up to 21 days later (p-value for standard ANOVA = 1.75 × 10 −5 ; omega squared = 0.193) (except for Tlen01 and Tatr01). The effect size for the shoot height after 7 days, as well as after 14 days and 21 days, was large; the observed differences among strains for the shoot height were of practical significance. At day 21 in the greenhouse, the lengths of shoots from seeds coated with Trichoderma strains were 4% to 11% longer, depending on the isolate, than those of untreated seeds. The best stimulatory effects on seedling growth were observed with the isolates Tahz02 and Tlix01 compared to the other tested strains. These effects were manifested by a highly significant (p-value for Kruskal-Wallis rank sum test < 0.0001) increase in both shoot (rank epsilon2 = 0.644) and root (rank epsilon2 = 0.535) length as well as fresh (shoot: rank epsilon2 = 0.591; root: rank epsilon2 = 0.621) and dry biomass (shoot: rank epsilon2 = 0.602; root: rank epsilon2 = 0.661) (Figure 3a-f). For all these variables, rank epsilon2 values showed a large effect size between strains; the magnitude of the differences in the mean scores among the strains seems to be large on these variables. significant (p-value for Kruskal-Wallis rank sum test <0.0001) increase in both shoot (rank epsilon2 = 0.644) and root (rank epsilon2 = 0.535) length as well as fresh (shoot: rank epsi-lon2 = 0.591; root: rank epsilon2 = 0.621) and dry biomass (shoot: rank epsilon2 = 0.602; root: rank epsilon2 = 0.661) (Figure 3a-f). For all these variables, rank epsilon2 values showed a large effect size between strains; the magnitude of the differences in the mean scores among the strains seems to be large on these variables. Additional ecophysiological traits correlated with the nitrogen status, such as leaf chlorophyll content and nitrogen balance index (NBI), as well as the secondary compounds related to plant defenses with the epidermal flavonols, were measured using Dualex, a non-destructive method commonly used for assessing ecophysiological traits in Additional ecophysiological traits correlated with the nitrogen status, such as leaf chlorophyll content and nitrogen balance index (NBI), as well as the secondary compounds related to plant defenses with the epidermal flavonols, were measured using Dualex, a non-destructive method commonly used for assessing ecophysiological traits in many plant species, including wheat [104]. Chlorophylls are key compounds for the photosynthetic process, and a great deal of research revealed the direct relationship between the number of chlorophylls and photosynthetic rates. As for the NBI, it is defined as the ratio of chlorophylls to epidermal flavonols (Chl/Flav) [104]. This proxy is an informative indicator for crop growth as the high and low nitrogen status can be quickly monitored by measuring NBI, providing accurate information for farmers to make timely N management decisions. It is generally assumed that a high NBI is an indicator of both high nitrogen content and high nitrogen use efficiency [105]. Results showed a remarkable impact on total chlorophyll content when seed coating with Trichoderma strains (p-value for Kruskal-Wallis rank sum test < 0.0001; rank epsilon2 = 0.411) (Figure 3g). All tested Trichoderma strains (except for Tlen01) significantly enhanced leaf pigments as demonstrated by the measurement of chlorophyll content: Trichoderma-treated plants exhibited 107% to 112% higher chlorophyll content which was positively correlated with the nitrogen content [106]. The improvement in the photosynthetic capacity of plants, modulated by various endophytic Trichoderma spp., would be attributed to an increase in the number of photosynthetic pigments (in addition to the expression of genes regulating the biosynthesis of chlorophyll and proteins integrated into the light-harvesting complex or the Calvin cycle [107]). These cell responses would be elicited by specific secondary metabolites and/or VOCs emitted by the fungus [108,109]. All isolates secreted VOCs (phenomena described in Section 3.4.2), and it is plausible that these emissions are involved in the observed increase in photosynthetic pigment content and in the ensuing promotion of plant growth.
Similarly, NBI, as well as leaf flavonol content, were significantly influenced by the strains (p-value for Kruskal-Wallis rank sum test < 0.0001), and large effect sizes were detected (rank epsilon2 = 0.333 for NBI, and 0.369 for flavonol). NBI has significantly increased between 108% to 118% under the effect of all Trichoderma strains except for the strain Tlen01, which has almost the same NBI as the control plants (Figure 3h). In contrast, leaf flavonol content decreased significantly for Tahz02 (−8%) and Tlix01 (−5%) and showed a similar trend for Tahz01 (−28%). Tatr01 and Tahz03 and Tlen01 showed the same flavonol content as the control plants ( Figure 3i). Flavonols (or flavonoids) are abundant and ubiquitous secondary compounds in plants. They are produced through the phenylpropanoid pathway, and their abundances can be regulated through de novo biosynthesis or rapid translocation and modification of existing compounds. Because of their remarkably versatile function, these molecules integrate many physiological processes, including plant growth (as well as disease resistance, which will be discussed in Section 3.4.1.). In this respect, their accumulation is a rapid and long-lasting physiological response to PGP in general and Trichoderma spp. in particular [110][111][112]. In our study, the slight decrease in (only for Tahz01, Tahz02 and Tlix01) or stable abundances of these flavonols for the other strains are not contradictory. Indeed, the chlorophyll content is positively correlated with the nitrogen content, while the content of epidermal flavonols is inversely correlated to nitrogen content [106]. Thus, the decreases in flavonoid accumulations observed with the strains exhibiting the strongest PGP potentials would reflect basal constitutive amounts of phenolic compounds in wheat leaves that can remain under the influences of external stimuli, a priori of root origin and induced under Trichoderma activity. In other words, these differential contents between strains would result from molecular and systemic signals, which would be strain dependent. In addition, our study is based on the duration of the interaction and not on short time intervals after the protagonists' interactions. In this respect, the abundances recorded 21 days after germination do not prejudge those that could have been recorded after a few hours or a few days of treatment. For example, a low level of flavonols correlated with a high NBI was also observed in maize seedlings from seeds inoculated with T. atroviride [106]. Similarly, we hypothesize here that the low flavonol abundances recorded for some strains echo the plant growth promotion induced by these same strains, knowing that under conditions of no stress, these two physiological processes interfere with each other. In any case, the content of epidermal flavonols recorded in wheat plants treated with the six Trichoderma isolates confirms that treatment with these strains did not affect the general health of the plants.

Trichoderma Isolates Exhibit Differential Biochemical Potentials
The isolated Trichoderma strains in this study developed significant and long-lasting PGP abilities four weeks after seed coating, and these abilities are strain-dependent. This finding strongly suggested that these strains produced growth-regulating factors or displayed metabolic activities that are involved in the increase of seedling growth kinetics. In this respect, the PGP potential of the six Trichoderma strains was assessed in vitro through the production of ammonia (NH 3 ), Hydrogen cyanide (HCN) and indole-related compounds and their ability to solubilize phosphate (Ca 3 (PO 4 ) 2 ) ( Table 2). Values are means ± SE (n = 3). Distinct letters in the same assay indicate significant differences (p ≤ 0.05). A "-" corresponds to an absence of phosphate solubilization and HCN or indole-like compound productions.
In our study, only the ammonia and the indole-related derivatives tests were positive. All strains were able to produce ammonia at the same level, and a level in the range of 14.1 µg/mL to 16.7 µg/mL. Tahz01 was found to produce the highest amount of ammonia. This result is supported by several studies showing the ability of Trichoderma to produce ammonia [61,113,114]. Ammonia is an essential trait linked to plant growth promotion, and the produced ammonia has been shown to supply nitrogen to the plant, thereby promoting root and shoot elongation and biomass [115]. Likewise, all strains produced indole-related derivatives, with the exception of T. lixii (Tlix01), for which no production was detected for the three technical repetitions. T. afroharzianum strains seem to be the most productive, in particular Tahz02 with 12.4 µg·mL −1 and, to a lesser extent, Tahz01 with 4.2 µg·mL −1 . The other isolates showed lower production ranges (i.e., 0.93 to 1.7 µg·mL −1 ). These levels were commonly recorded for some Trichoderma spp. [116][117][118] or rhizospheric PGPs [119,120]. Several factors can be put forward to explain these variations in the production of indole-like derivatives between strains: for example, the differential modulation of the biosynthetic pathways and the related-regulatory pathways or the thermodynamic specificities of enzymes that convert primary heterocyclic aromatic organic compounds (Tryptophan) into several indoles and related conjugates forms [121]. A colorimetric method using Salkowski's reagent was performed to quantify the production of indole-related compounds, but it does not allow the specific determination of Indole-3-acetic acid (IAA), even if the standard curves were carried out with pure IAA. Indeed, the levels of IAA produced in plant biological tissues or by microorganisms are infinitesimal to the order of nM [122]. However, Trichoderma spp. increased systemic biomass production and lateral root growth promotion in several plant species such as Arabidopsis [123], tomato [124] or cucumber [125], two physiological processes which are known to be dependent on IAA. We do not exclude that IAA is produced by our isolates, but this hypothesis must be validated by additional and more precise analytical methods such as High-Pressure Liquid Chromatography (HPLC) [126], Liquid Chromatography Electrospray Ionization Tandem Mass Spectrometric (LC-ESI-MS/MS) [127] or again by High-Performance Thin Layer Chromatography HPTLC [128]. It is nevertheless true that Salkowski's reagent test is a useful primary index in the screening of PGP functional rhizospheric microorganisms correlated to the production of various indole-related compounds and derivatives-among which IAA is a key substance-such as Trichoderma genus.
It is suggested that over 80% of rhizosphere bacteria are capable of synthesizing IAA and/or indole-related derivatives. Correlated to these specific metabolic syntheses, these strains are identified as biostimulators of plant growth, particularly by stimulating lateral and adventitious root length [129]. An increase in root surface area allows plants to access more minerals and nutrients from the soil and, presumably, to absorb them better [130]. It is confusing but nevertheless very interesting that T. lixii is one of the best-performing PGP isolates while it does not seem to produce indole-like derivatives. Although T. lixii is usually presented as the sexual state of T. harzianum in nature, we demonstrate here that each strain is genetically distinctive, including in their ability to produce Indole-like metabolites. To the best of our knowledge, no bibliographic reference reports the lack of production of indoles-like and/or IAA-related derivatives by T. lixii spp. Further analytical research (e.g., HPLC, LC, or HPTLC) would provide a definitive answer to this observation.
A contrario, all strains did not appear to metabolize phosphate nor produce HCN. Phosphate is essential for long-standing crop production, and although it is present in large amounts in soils, it needs to be mainly solubilized by the soil microbial community in order to be mobilized by plants [131] (Fankem et al., 2006). As for HCN, its production is intricately related to antifungal activity and the priming of the root length and root hair germination [132]. Each test was performed qualitatively on Pikovskaya's agar medium, which was supplemented with tri-calcium phosphate as an insoluble phosphate source, and by using Glycine or Succinate as primary substrates, respectively. The absence of these activities has been found in other studies applying the same methods [61,113,118], although some positive activities have also been demonstrated [133,134].

Trichoderma Isolates Exhibit Functional BCA Abilities against F. culmorum
Trichoderma spp. uses several antagonistic strategies against plant diseases. This includes the expression of a direct aggressiveness against its prey (i.e., the phytopathogenic fungi) by using mycoparasitism, antibiosis and secretion of extracellular cell wall-degrading enzymes and by eliciting transcriptional expression of a large set of host resistance-related genes.

Trichoderma Isolates Annihilate the Expression of Fusarium Seedling Blight Disease (FSB) Symptoms in Wheat
In order to confirm the bioprotective capacities of our isolates against Fusarium seedling blight (FSB) induced by F. culmorum in wheat, two bioassays were conducted, including a germination test in a water agar medium (in vitro), followed by a greenhouse seedling growth experiment for up to 3 weeks.
The first bioassay, based on seed emergence in vitro, showed highly pathogenic behavior of F. culmorum, with mycelium surrounding the seed by forming thick and robust hyphae ( Figure S4 in Supplementary Materials). This infection blocked seed germination before the appearance of coleoptile and radicle, killing more than 23% of the seeds (Figure 4a). Therefore, the final germination percentage (FGP) was significantly impacted by F. culmorum; it ranged from 75% in seeds infested with F. culmorum compared to 97% in non-infested seeds (control) (Figure 4a). For seeds that managed to germinate, they were significantly impacted by the aggressiveness of F. culmorum, showing the lowest total fresh weight and reduced root and shoot lengths (Figure 4a Treatment with Trichoderma isolates showed that they were able to profoundly combat the negative impacts of F. culmorum, which is reflected in the suppression of the appearance and development of any Fusarium symptoms ( Figures S3 and S4 in Supplementary Materials). This results in a significant increase in the final germination percentage for Tahz02, Tatr01, Tlix01 and Tlen01 to levels similar to that of the control and complete suppression of the negative effect of F. culmorum on total fresh seedling weight, shoot and root length for all strains. In addition, Tahz02, Tahz03, Tatr01 and Tlix01 showed a bioprotective and biostimulatory effect at the same time (in the presence of F. culmorum) by significantly increasing root length compared to the control (Figure 4b).
The second bioassay carried out in a greenhouse aimed to confirm the protective performance of Trichoderma isolates on growing wheat. Results showed that F. culmorum severely negatively affected the survival of seedlings (Figure 4d). Only 22% of seedlings had survived and showed reduced growth parameters (shoot/root lengths, fresh/dry total weights, leaf Chl/Flav content and NBI). We note, however, that the protective abilities of Trichoderma isolates were expressed differently between our six strains. The strain Tahz01 had completely protected the seedlings from F. culmorum, with neither the survival percentage nor the growth parameters being significantly influenced by F. culmorum infection, allowing seedlings to grow normally. In this respect, Tahz01 is considered to be the best protective Trichoderma strain tested in this study. In addition to that, this strain (Tahz01) has exhibited a biostimulatory effect in addition to its protective role. This effect was made very clear from the values of the shoot/root lengths, total fresh weight, shoot dry weight and leaf chlorophyll contents, compared to the values of the control seedlings (Figures 3 and 4d). The rest of the strains were partially protective by reducing the aggressiveness of F. culmorum. Tahz02 and Tahz03 were the most efficient when it came to seedling survival, while Tatr01, Tlen01 and Tlix01 were less effective (Figure 4d). Concerning the growth parameters, the ecophysiological traits recorded in leaves (chlorophyll content and NBI proxy) were relatively unaffected, while a weak protective action was noted in roots (fresh and dry weights) for Tahz03, Tatr01 and Tlen01 (Figure 3). infested seeds (control) (Figure 4a). For seeds that managed to germinate, they were significantly impacted by the aggressiveness of F. culmorum, showing the lowest total fresh weight and reduced root and shoot lengths (Figure 4a  Treatment with Trichoderma isolates showed that they were able to profoundly combat the negative impacts of F. culmorum, which is reflected in the suppression of the appearance and development of any Fusarium symptoms ( Figures S3 and S4 in Supplementary Materials). This results in a significant increase in the final germination percentage for Tahz02, Tatr01, Tlix01 and Tlen01 to levels similar to that of the control and complete Finally, plants that were diseased by inoculation of F. culmorum alone showed a drastic drop in the content of epidermal flavonols (Figure 3i). Interestingly, flavonol abundances were significantly up-regulated when the Trichoderma were co-inoculated with F. culmorum, except for Tatr01 and Tlen01, whose levels stayed stable and similar to the control. It is plausible that these increases in flavonol productions have a significant impact on the resistance of wheat induced by these Trichoderma, as it has been observed in other positive biotic interactions involving Trichoderma spp. or symbionts [110][111][112]135,136]. Finally, the synthesis of these secondary metabolites involves the intervention of several specific enzymes (phenylpropanoid pathway). Their regulation also includes accumulation levels, which are ensured by the transcriptomic regulation of their related encoding genes. These genomic regulations will be discussed in greater detail in the next item.

Trichoderma Display Antagonistic Abilities In Vitro against F. culmorum Growth
The antagonistic ability of the six native strains against F. culmorum was evaluated using in vitro dual culture assay. All the evaluated strains were capable of reducing F. culmorum growth. The inhibition rate (IR) was above 50%, indicating that all of the isolates exhibited high antagonistic activity. The most effective strain was Tlix01 with 73.8% IR, and that was only significant with Tahz03, which showed the lowest IR (57.8%) ( Figure 5; Figure S5 in Supplementary Materials). In addition, the results of the determination of diffusible metabolites demonstrated that the growth of F. culmorum was inhibited by at least 43.5%. This could be explained by the fast growing ability of Trichoderma to quickly colonize the surface and outcompete for resources, producing metabolites that lead to the suppression of the growth of F. culmorum [137].  The method used to determine volatile metabolites proved that all tested Trichoderma strains were able to inhibit F. culmorum development by 25.2% to 51.8%, with a significant difference between Tatr01 and Tahz03, respectively (Figure 5a,b; Figure S6 in Supplementary Materials). Trichoderma species produce various secondary metabolites, such as dif- The method used to determine volatile metabolites proved that all tested Trichoderma strains were able to inhibit F. culmorum development by 25.2% to 51.8%, with a significant difference between Tatr01 and Tahz03, respectively (Figure 5a,b; Figure S6 in Supplementary Materials). Trichoderma species produce various secondary metabolites, such as diffusible and volatile organic compounds (VOCs) [138], and among them, a group of more than 500 carbon-based compounds is able to diffuse into the atmosphere and soil and/or is solubilize in water. These compounds are known as microbial volatile organic compounds (mVOCs), and their biosynthesis is dependent on strains and environmental conditions [24,108]. They display multiple biological functions, such as enhancing plant growth and resistance [108,139,140], and are harmful to a broad spectrum of plant pathogens [24,29,140] or both [141]. The molecular diversity of VOCs is remarkable; it includes (thio)alcohols, aldehydes, heterocycles, hydrocarbons, ketones, phenols, thioesters and a plethora of derivatives [142,143]. Our study has demonstrated the ability of the Trichoderma isolates to produce volatile and soluble mVOCs, which would weaken F. culmorum mycelia development to be vulnerable to hydrolytic enzymes. These mycoparasitic events would lead to the successful inhibition of F. culmorum growth, as previously observed in the interaction Trichoderma-Pyrenophora teres, the causal agent of barley net blotch [144]. In addition, we do not exclude the possibility that some of the emitted mVOCs are involved in the biostimulation of seed germination and seedling growth, especially in the presence of F. culmorum, as the general metabolism of Trichoderma is influenced by the roots of its host [145,146], as well as its potential prey including Fusarium graminearum [147]. This hypothesis would partly explain the PGP ability of Tlix01, where mVOCs could compensate for the absence of IAA-related indolic metabolites in this strain.
Along the same lines of exploring Trichoderma weapons used in the control of F. culmorum development, qualitative determination of enzymes, presumably involved in the antagonist activity of these isolates, revealed positive results illustrated by bubbling for Catalases and a halo of degradation inside the mycelial growth zone for Endoglucanases, Proteases, Chitinases and Amylases production ( Figure 6).
Microorganisms 2023, 11, x FOR PEER REVIEW 20 of 36 [147]. This hypothesis would partly explain the PGP ability of Tlix01, where mVOCs could compensate for the absence of IAA-related indolic metabolites in this strain. Along the same lines of exploring Trichoderma weapons used in the control of F. culmorum development, qualitative determination of enzymes, presumably involved in the antagonist activity of these isolates, revealed positive results illustrated by bubbling for Catalases and a halo of degradation inside the mycelial growth zone for Endoglucanases, Proteases, Chitinases and Amylases production ( Figure 6). These enzymes are produced as part of Trichoderma metabolism to obtain nutrients and to defend themselves against other microorganisms. Catalases, which are oxidative enzymes, have been shown to be involved in its biocontrol activity against F. culmorum, while the hydrolytic enzymes are secreted to break down Fusarium cell walls to get access to nutrients and help thus to reduce the severity of Fusarium infections. We hypothesize These enzymes are produced as part of Trichoderma metabolism to obtain nutrients and to defend themselves against other microorganisms. Catalases, which are oxidative enzymes, have been shown to be involved in its biocontrol activity against F. culmorum, while the hydrolytic enzymes are secreted to break down Fusarium cell walls to get access to nutrients and help thus to reduce the severity of Fusarium infections. We hypothesize that these lytic enzymes (which are primarily involved in the extracellular digestion of Trichoderma) act in concert with various antagonistic metabolites (mVOCs) to take part in the mechanisms that drive Trichoderma's hyperparasitism activities against F. culmorum.

Trichoderma Strains Stimulate the Wheat Immune Responses General Consideration
As for the biostimulation of plant defenses, all isolates significantly modulated the transcriptional expressions of genes involved in the system-wide immunity of wheat (all p-values < 0.05) (Figure 7; Table S2 in Supplementary Materials). Plant defenses are remarkably diverse, and the ways in which they are regulated are also remarkably varied and complex. They involve multiple pathways (viz., PR proteins, secondary metabolites, parietal reinforcements, reactive oxygen species (ROS)) inherent to the recognition of bioagents, which, in our study, correspond to both pathogenic (F. culmorum) and beneficial agents (Trichoderma spp.). The accumulation of plant defenses is one among cellular responses that are set off by the recognition of external agents or substances and the self-recognition of plant molecules released by itself during the interaction such as stress phytohormones (SA, JA, Et, IAA, etc.), ROS with the hydrogen peroxide (H 2 O 2 ), and/or various metabolites coming from microorganisms living on the host plant. The defense markers we selected for this study are common and cover several defense pathways (PR-proteins, phenylpropanoids, sesquiterpenoids, oxidative stress, recognition and signal transduction).
Molecular analysis was carried out on the roots and leaves of 21-day-old seedlings grown from seeds that were coated with the six Trichoderma isolates. An overview of transcriptional data showed that all targeted molecular markers were significantly modulated by the application of fungus (Trichoderma/Fusarium) independent of analyzed tissues. The modulations are highly diversified, with positive and negative kinetics reflecting the systemic regulation of the defense responses established by the host, which perpetuate over time during plant growth. These physiological features highlight the complexity and richness of the molecular responses that the plants deploy to provide a general resistance phenotype correlated to the BCA capacities of Trichoderma or the disease susceptibility caused by the virulence of F. culmorum.
Another general observation is that Trichoderma, when co-inoculated with F. culmorum, significantly enhanced the expression of some of the positive modulations that were recorded in the case of inoculation with only Trichoderma or even completely reversed other transcriptional expression trends induced by F. culmorum when inoculated alone (downvs. up-regulations, and vice versa). These cellular responses clearly demonstrated that Trichoderma strains triggered deep and long-lasting transcriptional reprogramming of genes involved in plant immunity of wheat seedlings that were infected with F. culmorum. The reprogramming concerned all cellular processes: recognition (e.g., LRR, RLK, RLCK and WAK6), transduction pathways (CERK3), oxidative burst (ROS generation with RbohD and -F, and the antioxidant system with GPX and CAT), accumulation of defenses (secondary metabolisms with PAL, CHS, GST and COMT, as well as PR-proteins with DEF and LTP), and modulation of phytohormonal balances (ACO and AOS2).
recognition of plant molecules released by itself during the interaction such as stress phytohormones (SA, JA, Et, IAA, etc.), ROS with the hydrogen peroxide (H2O2), and/or various metabolites coming from microorganisms living on the host plant. The defense markers we selected for this study are common and cover several defense pathways (PR-proteins, phenylpropanoids, sesquiterpenoids, oxidative stress, recognition and signal transduction). Defense-related gene modulation in wheat roots and leaves following fungal inoculations with Trichoderma isolates alone or in combination with F. culmorum. Organs sampled for molecular analyses came from plants used for physiological and morphological analyses. Transcript levels for each gene were estimated using real-time qRT-PCR analyses and normalized by the expression of three housekeeping genes. Relative transcript abundance rates were obtained by the 2 −∆∆CT method. Data correspond to the mean of three independent biological experiments; each biological repetition corresponds to a mix of 5 plants. Colors represent the levels of transcriptional expressions (blue, up-expression; red, down-expression, compared with healthy samples treated with sterile PDB). NOX, NADPH oxidase; AQP, Aquaporin. Statistical analyses were detailed in Table S2 in Supplementary Materials. Gene acronyms and related primers used for qPCR amplification are given in Table S1 in Supplementary Materials. SA, salicylic acid; Et, Ethylen; JA, Jasmonates.
When analyzing transcriptional responses at the Trichoderma isolates level, it is interesting to note that Tatr01 and Tlen01 exhibited contrasting profiles to those observed with the rest of the isolates. This was significant for genes involved in the recognition process and transduction pathways (e.g., WAKs, MAPK3, WSF and NFLX), oxidative stress (RbohF, AQP, GPX, SOD and GST), defense metabolites (PAL, FLS, PPO, CCR3, GSL22, SQTS, CHS, CCR3, GSL22, PPO, PR1b, GLU, and OXO) and phytohormonal balances (mainly Etand JA-dependent pathways with ACO, LOX, AOS1 and -2). These two isolates (Tatr01 and Tlen01) were protective, similar to the other isolates, against the emergence of morbid disease symptoms (root and leaf rots) caused by F. culmorum (Figure 4b; Figure S3 in Supplementary Materials) but not against growth delays due to the presence of F. culmorum on plants (Figure 3). Although these are correlations, this finding may be explained by the fact that Tatr01 and Tlen01 have affected the transduction pathways inherent to plant growth and immunity in a different way than other isolates. This data is a prerequisite for guiding the selection of strains to be developed in the field. To our knowledge, such contrasting strain-dependent modulations in the instruction of different physiological functions are rarely demonstrated or further discussed in wheat interacting with BCAs in general and with Trichoderma in particular.

Trichoderma Isolates Stimulate the Systemic Accumulation of Wheat Defenses
Plant defenses are expressed sequentially, temporally and spatially. Trichoderma inoculation in the presence or absence of F. culmorum modulated all the targeted markers. Based on our transcriptomic results, we proposed a putative model of the long-term molecular interaction between above-ground and underground parts of wheat seedlings grown from seeds coated with the six Trichoderma isolates.
The expression levels of the defenses encoding genes resulting from the recognition of the agents by plant receptors (R) would participate in triggering both the root colonization of Trichoderma and the systemic modulation of defenses accumulation which then leads to the resistance of wheat against F. culmorum. Different molecular modalities of Trichoderma spp. plant recognition and its related defense elicitation have been described, for example, involving the Sm1/Epl1 elicitors [148,149]. Consistent with several transcriptomic and proteomic studies [30,150], the six Trichoderma isolates induced the modulation of several plant resistance (R) genes, which was correlated with the activation of a wide range of plant defenses and a high physiological level of resistance. It is then also plausible that this resistance acquisition should be extended to pathosystems other than F. culmorum (and of which WAK2 would be one of their recognition systems [151]), such as Blumeria graminis f. sp. tritici (LRR RLK) [152], Puccinia triticina (WAK6) [153], Rhizoctonia cerealis (RLCK) [154], and those whose related receptor-like kinases encoding genes are significantly up-regulated by certain Trichoderma strains alone and/or in the presence of F. culmorum. However, the enhancement of the effector-triggered resistance relies on the highly specific interaction between the microorganism's effectors and the potential plant's receptors (R). The fact that we observe diversified transcriptional patterns of the (R)-gene responses between strains suggests that each isolate has deployed specific effectors that interfere with the plant effector-triggered immunity (ETI). It will be very informative to further characterize the effectors involved in the interaction of the most efficient Trichoderma strains on wheat.
Receptors perceive and process signals from invading microorganisms in diverse pathological and beneficial systems. Among the subsequent cellular responses induced by our Trichoderma isolates, we observed an increase in transcription of wheat respiratory burst oxidase homologues (i.e., NADPH oxidase with the Rboh-D and -F subunits) and proteins involved in ROS diffusion (AQP) and scavenging (antioxidative complex). These modulations mainly involved up-regulations that led to the modification of the global oxidative state of plants that were linked to a generation of ROS. They occurred in the roots (except for Rboh-D) and in the leaves, where they were correlated to the modulation of the SA-dependent markers. Indeed, SA is a key mediator in eliciting an oxidative burst in challenged plants, thus regulating antioxidant metabolism. In this respect, these results agree with those which were reported in wheat inoculated with different Trichoderma species where the application of Trichoderma triggers both oxidative stress and systemic defense pathways in wheat seedlings [155].
Besides the toxic and antimicrobial aspects of ROS, hydrogen peroxide (H 2 O 2 ) is a key player in the complex signaling network of plant responses to abiotic and biotic stresses. It is, therefore, essential to maintain ROS at sub-toxic levels, which reflects a delicate balance between ROS production (involving the activation of ROS-generating enzymes under stress and the unavoidable production of ROS during basic cellular metabolism) and ROS elimination pathways. Here, our experiments demonstrated that Trichoderma induced persistent and latent oxidative stress in wheat, but this stress was counteracted by the production of antioxidant enzymes (CAT and SOD) that then reduced the oxidative damage linked to the high chemical reactivity of ROS. In addition to their role in the elimination of F. culmorum, these ROS and in particular H 2 O 2 , would also participate in the changes of the plant's phytohormone balances rather than in the elimination of F. culmorum, which, subsequently, would notably activate the accumulation of defenses. Furthermore, H 2 O 2 is a relatively stable molecule over time, and it has the distinctive feature of freely diffusing through the plasma membrane. However, its diffusion in tissues should and could be regulated by peroxyporins. Peroxyporins are specific aquaporins (AQP) that belong to the Major Intrinsic Protein (MIP) family [156]. Aquaporins are studied in a plethora of cellular physiological pathways due to their ability to regulate water transport in biological systems. However, peroxyporins play a relevant role in controlling H 2 O 2 permeability, and they ensure the ROS outcome during oxidative stress, including that related to the initiation of the disease immune responses by plant cells [157]. AQP expressions were modulated by all isolates, both in roots and leaves, implicating them in the diffusion of many solutes and also H 2 O 2 , reinforcing the hypothesis that ROS-dependent signals may have influenced the overall resistance responses of the plant. It is, therefore, interesting to correlate the cellular pathway inherent to the production, diffusion and pro-oxidative control of ROS with the resistance of wheat against F. culmorum by Trichoderma.
Concomitant with oxidative stress, the strong induction of various defense pathways is a key plant defense strategy to limit pathogen development. It involves the production of a large panel of "biopesticide" secondary metabolites, including phenylpropanoids, terpenoids (SQTS, FPS1), and enzymes involved in the biosynthesis of monolignol (Cinnamyl alcohol dehydrogenase, CAD, and Cinnamyl coenzyme reductase 3, CCR3), which participate in cell wall strengthening and callose apposition (GSL22). These genes are mainly induced in leaves and roots, and a few of them may also be repressed but with similar modulation rates.
PAL and CHS are key enzymes of the phenylpropanoid biosynthesis pathway, providing precursors for the synthesis of antimicrobial phenolics such as lignins, flavonoids, isoflavonoids and coumarins. Furthermore, these enzymes are de facto involved in the SA defense pathway since SA is synthetized from phenolic derivatives. However, more broadly, it should be emphasized that several genes involved in the biosynthesis of flavonoids and terpenoids were up-regulated by several Trichoderma isolates, thus reinforcing the relevance of such a pathway during the Wheat-Trichoderma interaction leading to efficient systemic resistance against F. culmorum. These transcriptional up-regulations are correlated with the leaf epidermal flavonol contents recorded in the leaves of plants treated with both Trichoderma isolates alone or with F. culmorum (Figure 3i).
The intertwining of enzymatic markers involved in different metabolic pathways under stress seems to contribute to disease resistance. For instance, CAD and CCR3, two key enzymes involved in lignin biosynthesis, interplay in apoplasm with the oxidative enzymes Peroxidase (POX) and Polyphenol oxidase (PPO) in the presence of ROS (O 2 •− , • OH) and H 2 O 2 to catalyze the formation of lignins and their crosslinking with other oxidative phenols, proteins and molecular components of the cell wall. These biochemical interactions contribute to the reinforcement of the cell wall and to the callose apposition process itself, assisted by Glucan synthase-like activities (GSL22). Interestingly, these key enzymes display significant up-coregulation in the roots and leaves elicited by Trichoderma alone or with F. culmorum. In addition, their accumulation patterns are analogous to Superoxide dismutase (SOD) enzymes that are responsible for the dismutation of superoxide radicals into oxygen and H 2 O 2 . A contrario, the expression profile of these enzymatic markers was significantly inhibited in the roots of Trichoderma-F. culmorum co-inoculated plants, compared to F. culmorum alone. These expressions are possibly related to the toxic impact of certain enzymes on the development of microorganisms in general and on the longterm Trichoderma colonization in roots in particular. For example, PPO activity oxidizes phenols into quinones, which are more toxic to microorganisms than the original phenolic compounds [158]. Likewise, POXs belong to a large multigene family that participates in several physiological processes, such as the synthesis of phytoalexins, lignin and suberin, the cross-linking of cell wall components and the rapid synthesis of ROS derivatives by oxidative burst contributing to cell death [159]. All these molecular responses provide a physical barrier to and/or substantially limit pathogen spread and invasion throughout the plant. Finally, their inhibition in co-treated root correlates with the down-regulation of the SA-and JA-dependent signaling pathways, possibly to finely control the abundance of related defenses in the roots which are in direct contact with the fungal protagonists.
As for the pathogenesis-related (PR) proteins, which are grouped into 17 families based on their sequence and biological activities, they are components of an early response toward pathogen infection, playing a key role in pathogen-spread limitation [160]. Phytohormones (SA, JA and Et) with ROS are involved in the signal transduction network and the accumulation of PR-proteins primed by Trichoderma spp. [32,[161][162][163]. In our assay, PR-protein transcriptional accumulations were significantly activated in leaves challenged by Trichoderma alone, while these expressions were strongly potentiated in the leaves and also in roots in the case of co-treatment Trichoderma-F. culmorum. Finally, co-regulations of PR-proteins, as well as other defense markers, could be pointed out mainly in the leaves challenged by Trichoderma alone; it concerns CHTb, POXa, DEF, LTP, OXO and GLP, with LRR-RLK, CERK3 and GPX. Interestingly, these expression patterns are quite connected to the modulation of the JA-dependent signaling pathways (AOS1 and AOS2), knowing that among these markers CHTb (PR-4), GLP and AOS2 are also positively induced by fungal pathogens through the JA signaling pathway [164]. To our knowledge, these markers do not seem to be linked to a common signaling pathway. However, the fact that these profiles are correlated may underline key regulatory points in the establishment of certain disease resistance, such as those involved in the role of the LRR-RLK receptor against Blumeria graminis f. sp. Tritici, the agent of wheat powdery mildew.
Finally, we do not omit that some defenses are weakly/unmodulated or significantly inhibited during the interaction in roots in similar patterns to those recorded with F. culmorum alone (e.g., C4L, CCR3, GSL22). It is plausible that these defenses remain under the control of the pathogen, even with the strong regulations that Trichoderma exerted on the global signaling pathways in wheat seedlings. However, these modulations could also be related to the successful, long-lasting Trichoderma establishment on or in these plant organs. We can then assume that certain specific defenses are "effective" against F. culmorum (correlated to their down-regulation) and that they would be equally toxic to Trichoderma. For example, the expression profiles of certain PR-proteins (PR1b, GLP, OXO), key biosynthetic steps of secondary metabolites or the sextet (SOD-PPO-POX-CCR3-CAD-GSL22) would fit this hypothetical scenario. As demonstrated, even though these defenses can be upregulated rapidly by Trichoderma treatments [30,116], their expression is reprogrammed over time to perpetuate an optimal mutualistic interaction that the plant maintains with its beneficial microorganisms. From this perspective, the mechanisms deployed by Trichoderma to persist in the root system as an avirulent symbiont are not yet fully revealed. It is highly probable that it depends on "virulence effectors" involved in plant defense suppression, whose nature and mode of action remain to be deciphered.

Trichoderma Isolates Influence the Wheat Phytohormone Balances
Root penetration by Trichoderma is reported to transiently induce the expression of defense-related genes and the production of antimicrobial compounds. Our results demonstrated here imply that root colonization by Trichoderma after seed coating application maintained the transcription of many defense-related genes at significant levels for a relatively long period of time (21 days). These observations suggest that the long-term response of wheat to Trichoderma may involve phytohormone signaling.
In addition to the primary mechanisms of elicitor recognition by plant receptors, complex phytohormonal adjustments occur that supplement or potentiate the primary signals of direct microbial recognition. It is unanimously recognized that beneficial agents in general, and Trichoderma in particular, have the potential to profoundly influence the phytohormonal balances of plants involved in local and systemic signaling responses to stress. These balances include salicylic acid (SA)-and jasmonate (JA)/ethylene (Et)-mediated signaling responsible for the Systemic Acquired Resistance (SAR) and Induced Systemic Resistance (ISR), respectively, making plants resistant to a broad spectrum of pathogens [116,165,166]. Generally, SAR is induced by fungal biotrophic pathogens, and ISR usually relies on the interaction of plants with symbionts or necrotrophic microorganisms. However, Trichoderma spp. are reported to engage several hormonal-signaling pathways to instruct the elicitation of disease resistance in plants. They act in synergy and in balance with each other, and they depend closely on the strains, the plant genotypes, the environment (including the challenged pathogen) and the temporal modalities of the interaction [116,167,168]. We explore these balances here through the SA-and JA/Et-mediated signaling pathways.
Concerning the SA-dependent pathway, Trichoderma isolates the negatively modulated NPR1 and ICS in the roots, whether inoculated alone or with F. culmorum. These modulations are quite diversified in the leaves, being positive for simple Trichoderma inoculation or with inhibitions or activations for co-inoculation with F. culmorum, but they are straindependent. The down-regulation of the SA pathway markers recorded in our work should not be opposed to the key role played by SA in the early stages of root colonization by beneficial microorganisms where SA accumulation and/or transcriptional up-regulation of SA pathway-dependent genes are reported. Early phases involve plant responses that resemble those induced by biotrophic pathogens present in the intercellular spaces, similar to many beneficial agents [169]. These first phases are correlated with transitory accumulations of SA that may be responsible for the accumulation of plant defenses. In an interaction involving Trichoderma, these common physio-molecular responses would prevent the fungus from entering the apoplastic space of the epidermal and cortical tissues towards the root vascular system [170]. Here, our study targeted the defense responses of 21-day-old wheat seedlings. This specific temporality of our analyses would then explain the weak modulation and inhibition of several defenses observed in roots and could also be related to the successful and long-lasting establishment of Trichoderma on these plant organs. Such partial suppression of SA-dependent responses in plants also appears to be necessary for the onset of symbiotic associations [171,172]. Moreover, these organ-and strain-specific transcriptional responses are very informative. They highlight that isolates behave and/or are perceived differentially by the plant, which is dependent on their presence both alone or with their plausible prey (F. culmorum). This demonstrates the complexity of the interactions that beneficial agents establish with plants in the presence or absence of other protagonists and whose outcomes are closely dependent on the biotic and abiotic environments where these interactions take place.
The expression pattern of the JA/Et pathway markers allowed us to discover some interesting correlations. Ethylene (Et) is a phytohormone that regulates plant growth, development and senescence, but low Et concentrations in the root system are correlated with higher shoot growth. These molecular responses (along with the ecophysiological traits we recorded) have to be linked to the fact that many PGP microorganisms from the rhizosphere could reduce the levels of 1-aminocyclopropane-1-carboxylic acid (ACC), which is the precursor metabolite of Et, by producing the enzyme ACC deaminase (ACCD). Trichoderma strains own this feature of Et regulation [173], regulating the level of ET in conjunction with other phytohormones involved in the growth process (GA3, CKs) [124,174]. This could be explained by the fact that Tatr01 and Tlen01 have induced an overexpression of ACO genes, which is correlated with a reduction in root growth (Figure 3d-f). In our analysis, it is interesting to note that the expression of ACS and ACO genes, which encode for enzymes involved in the last two key steps of the Et biosynthesis in plants, is apparently desynchronized between biological variations (i.e., inoculation alone or combined with F. culmorum). Similarly, the expression of regulatory factors known to be related to ethylene (Et)-mediated signaling (ERF3, PIE) have shown highly divergent profiles. These contrasting transcriptional profiles have not been over-interpreted to avoid making wrong conclusions, but they could plausibly be attributed to several global regulations. These regulations are generated from long-term transcriptional reprogramming of Et-signaling influenced by the general phytohormonal balances of plants, to which Trichoderma contributes with its own phytohormone-like production.
The same goes for the jasmonates (JA)-mediated signaling pathway. Jasmonates are a family of molecules derived from the lipid metabolism and subsequent peroxidation of membrane lipids caused by an oxidative burst. Jasmonates, in cooperation with other phytohormones (including SA), play an active role in the regulation of plant innate immune responses by controlling a cascade of complex physio-biochemical changes (metabolite accumulation and enzymatic activities) within the host plant, which then contributes to the development of resistance/tolerance against a wide range of invaders [175,176]. In addition, jasmonates influence the regulation of mutualistic interactions between plant and their associated microbes in the roots and aerial parts. In this respect, the fact that Trichoderma would inhibit JA-dependent regulatory pathways (LOX, AOS1 and -2) in the roots of wheat plants further mirrors the fact that mycorrhizal colonization is stimulated by low concentrations of jasmonates rather than higher ones. Concerning leaves, plant systemic resistance induced by Trichoderma spp. correlate with a significant increase in JA-related markers, which is consistent with several works [116,177]. However, a complete re-organization of the JA-marker transcripts accumulations occurs when Trichoderma is coinoculated with F. culmorum, while F. culmorum alone displayed inhibitor activity on these JA-markers. These trends suggest the possible existence of cooperative interactions between the mutualistic plant symbiont (i.e., Trichoderma) and the roots of the host plant in response to F. culmorum infection through the emission of specific JA derivative metabolites such as the methyl jasmonates (MeJA). Both molecules are natural inducers of plant resistance and their concomitant application with Trichoderma spp. significantly improves the induction of plant immune responses [175,177]. In addition, MeJA also exhibits significant antimicrobial activities by inhibiting spore germination and mycelial development of various plant pathogenic fungi [178][179][180].
These modulations need to be analyzed in a more holistic sense. We have to mention that most of the studied Trichoderma species colonize either the root surface or the interior as endophytes [181,182], including in wheat [183]. However, within the soil fungal communities, Trichoderma spp. has a relatively low abundance in both agricultural [184,185] and natural ecosystems [186,187]. The act of seed coating guarantees the presence of Trichoderma in wheat tissues over time to ensure a long-lasting and systemic elicitation of plant responses. Among the early signaling-induced events related to the host-Trichoderma recognition, high levels of SA, Et and JA were observed in different parts of the plant [188,189].
While several studies have focused on the early (i.e., hours or days) and systemic elicitation of Trichoderma when inoculated on roots and/or leaves by transcriptomic or proteomic analysis, our study has the originality of providing new insights about these root and leaf responses when seeds are coated with Trichoderma. The isolates thus have time to acclimatize and settle in their edaphic environment while potentially colonizing the growing plant organs (i.e., coleoptiles, then all or part of the young seedling, including the roots). The result after 21 days of interaction in a greenhouse shows that plant response patterns related to phytohormonal defense regulations are highly diversified between pathways (SA and JA/Et) and closely dependent on the applied strains. These responses are very complex to interpret, but growth and immune performances appear to be associated with the simultaneous transcript expression of SA and JA/Et defense-related genes, as well as with the expression of an oxidative burst and accumulation of antimicrobial metabolites. SA-and JA/Et-dependent signaling are often depicted as antagonistic, but here, their synergistic overlap seems to occur instead, as previously reported for some pathosystems [116,190,191]. However, although the induced SAR and IRS seem compatible in our study, the signaling pathways leading to plant resistance and growth are different and distinctly strain dependent. Furthermore, the complexity of the transcriptional responses observed would plausibly be linked to the interaction of two pools of phytohormones emitted during plant/Trichoderma interaction, i.e., the first emitted by the plant in response to the recognition and settlement of Trichoderma in the root, and the second emitted by the fungus to manipulate its (biotic) environment, and which, here, would fully integrate the signaling pathways of the plant. In other words, the global and systemic hormonal regulations observed in this study would be conditionally linked to the ability of Trichoderma to produce several phytohormones (such as ABA, CKs, GAs, SA, Et and IAA). The level of these phytohormones is strain-specific, and their cooperation synergistically modulates plant disease resistance and growth promotion [192].
Therefore, and consistent with what has been observed in other studies on tomato [193], cucumber [194] or olive [32], this evidence further underlines the key influence of Trichoderma on the overall phytohormonal balances that interconnect plant defense and development responses. The transcriptional expression levels of the markers studied here did not allow us to identify which phytohormone-likes are mainly produced by the different Trichoderma strains, nor the way in which they cooperate to reach this effect on plant growth and defense promotion. Other ongoing research will quantify the phytohormone-like compounds produced by these beneficial isolates in order to identify their potential mode of action during the molecular dialogues of the defensive protection process in wheat.

Long-Lasting Root Installation of Trichoderma Isolates for a Sustainable Plant Protection
It should be noted that the physiological responses recorded here are expressed in the roots and leaves of wheat plants harvested after three weeks under greenhouse conditions and that the fungal treatments were carried out only on seeds.
Several events can explain the durability of the signal over time. Each Trichoderma isolate was detected in the roots of plants coming from seeds treated with that same strain ( Figure S7 in Supplementary Materials). The same goes for the pathogen, which was only detected in the roots of plants originating or not from seeds previously treated with Trichoderma isolates. Regarding leaves, no PCR amplification signal was detected for all Trichoderma isolates and F. culmorum was detected in plants treated only with this pathogen. However, it is plausible that the isolates are also present on/in the aerial parts but at levels not measurable by the analytical approaches used here and/or that endophytic populations mask the signal. In this respect, these potential trace populations that would be controlled by the immune performance of plants stimulated by Trichoderma could conceivably still have a positive influence on the foliar defense responses by intensifying the primary systemic signals emitted by the roots during the interaction with Trichoderma isolates.
The prolonged presence of the fungal protagonists should thus have an influence (individually or synergistically) on the intrinsic "immune memory" of plant cells whose inhibitory feedback regulations have not yet been fully deployed. The combination of these two events can temporarily reduce the vitality of the growth of the plants. Plants actively prioritize defense over growth, and this prioritization and the associated growth trade-offs are regulated by crosstalk between plant hormones [195,196]. In our experimental context, these negative collateral effects were not observed with any of the beneficial interactions, and Tahz01, -02 and -03, and Tlix01 significantly improved plant growth dynamic ( Figure 3). Therefore, this data is crucial to determining the choice of strains whose positive effect is closely linked to their ability to develop BCA and PGP capabilities simultaneously.
Finally, in our study, where in planta experiments were conducted in a growth chamber, F. culmorum was traceable in diseased plants, as well as in symptomless plants that were previously treated with each Trichoderma isolate. This data is of great significance because the fact that Trichoderma persistently colonizes the roots alongside the pathogen suggests the relevance of using these agents as a form of biocontrol, particularly for a pathogenic population, and/or in the modulation of the expression of pathogenicity genes that instruct the level of pathogen's aggressiveness. The molecular dialogues that regulate cellular responses in each of the interacting fungal protagonists open up many remarkable exploratory avenues, which will be studied in the future.

Conclusions
To conclude, this study was conducted to isolate local Trichoderma spp. from rhizospheric soil of an old wheat variety (cultivated in an organic farming system) in order to use them as Plant Growth Promoting Fungi (PGPF) and Biocontrol Agents (BCAs) against F. culmorum, the agent responsible for Fusarium seedling blight disease in wheat (FSB). Of all the Trichoderma isolates, six strains were tested in this study. The results showed that seed coating with Trichoderma strains before sowing was efficient in promoting plant growth by increasing shoot/root length, total fresh/dry weight, total chlorophyll content and nitrogen plant status (NBI). In addition, strains appeared to be efficient biological control agents in their capacity to inhibit F. culmorum growth in vitro and by reducing disease severity under controlled environmental conditions. Moreover, the strains were able to induce systemic and long-lasting immune performance in the roots and leaves of 21-day-old seedlings. These physiological characteristics correlated with profound transcriptional changes in the gene encoding for plant defenses and reprogramming of hormonal balances, leading to simultaneous expressions of SA-and JA/Et-dependent pathways. The spatio-temporal coordination of these cellular responses conferred systemic resistance to wheat seedlings where the levels of its protective effect depended on the Trichoderma strain tested.
This research allowed us to explore the local Trichoderma strains which colonize the wheat rhizosphere and examine their performance as PGPF and BCAs. Their successful use here needs to be followed up by further analysis to fully understand the nature of the molecules (diffusible and volatile metabolites) involved in the biocontrol activity, as well as their impact on wheat growth and resistance performances in the field. In this respect, new studies focused on deciphering the molecular dialogue between the protagonists (Fusarium, Trichoderma and wheat) are in progress to better understand and valorize the isolated Trichoderma strains, which will then potentially guide their industrial application.
Supplementary Materials: The following supporting information can be downloaded at: https: //www.mdpi.com/article/10.3390/microorganisms11061512/s1, Figure S1: Tunisian site of origin of the Trichoderma strains used in this work; Figure S2: Macroscopic and microscopic aspects of the Trichoderma isolates; Figure S3: Root and leaves disease symptoms in wheat (Khiar variety) caused by Fusarium culmorum; Figure S4: Bioprotection effect of the six Trichoderma isolates against F. culmorum infestation on wheat ("Khiar" variety) in axenic environment; Figure S5: Direct mycelial confrontation between Trichoderma isolates and F. culmorum; Figure S6: Indirect antagonistic activity related to emission of volatile and diffusible organic compounds by Trichoderma isolates against F. culmorum; Figure S7: Qualitative molecular characterization of the presence of Trichoderma or F. culmorum isolates in the aerial and underground parts of wheat plants 21 days after seed treatment with Trichoderma and/or F. culmorum isolates; Table S1: Primers used for qPCR amplification of the wheat defenseencoding genes; Table S2: Statistical analyses of the transcript accumulation induced by Trichoderma isolates and/or F. culmorum.